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Abstract 

This paper considers the propagation of shallow-water solitary and nonlinear peri- 
odic waves over a gradual slope with bottom friction in the framework of a variable- 
coefficient Korteweg-de Vries equation. We use the Whitham averaging method, using 
a recent development of this theory for perturbed integrable equations. This general 
approach enables us not only to improve known results on the adiabatic evolution of 
isolated solitary waves and periodic wave trains in the presence of variable topography 
and bottom friction, modeled by the Chezy law, but also importantly, to study the 
effects of these factors on the propagation of undular bores, which are essentially un- 
steady in the system under consideration. In particular, it is shown that the combined 
action of variable topography and bottom friction generally imposes certain global 
restrictions on the undular bore propagation so that the evolution of the leading soli- 
tary wave can be substantially different from that of an isolated solitary wave with 
the same initial amplitude. This non-local effect is due to nonlinear wave interactions 
within the undular bore and can lead to an additional solitary wave amplitude growth, 
which cannot be predicted in the framework of the traditional adiabatic approach to 
the propagation of solitary waves in slowly varying media. 

1 Introduction 

There have been many studies of the propagation of water waves over a slope, sometimes 
also subject to the effects of bottom friction. Many of these works have considered linear 
waves, or have been numerical simulations in the framework of various nonlinear long-wave 
model equations. Our interest here is in the propagation of weakly nonlinear long water 
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waves over a slope, simultaneously subject to bottom friction, a combination apparently 
first considered by Miles (1983a,b) albeit for the special case of a single solitary wave, or a 
periodic wavetrain. An appropriate model equation for this scenario is the variable-coefficient 
perturbed Korteweg-de Vries (KdV) equation (see Grimshaw 1981, Johnson 1973a,b), 

A t + cA x + ^A+ %AA X + °^-A xxx = -C D -^\A\A. (1) 
2 2n o n z 

Here A(x, t) is the free surface elevation above the undisturbed depth h(x) and c(x) = 
a/ gh(x) is the linear long wave phase speed. The bottom friction term on the right-hand 
side is represented by the Chezy law, modelling a turbulent boundary layer. Here Cd is a 
non-dimensional drag coefficient, often assumed to have a value around 0.01 (Miles 1983a,b). 
Other forms of friction could be used (see, for instance Grimshaw et al 2003) but the Chezy 
law seems to be the most appropriate for water waves in a shallow depth. In ([1]) the first 
two terms on the left-hand side are the dominant terms, and by themselves describe the 
propagation of a linear long wave with speed c. The remaining terms on the left-hand 
side represent, respectively, the effect of varying depth, weakly nonlinear effects and weak 
linear dispersion. The equation is derived using the usual KdV balance in which the linear 
dispersion, represented by d 2 /dx 2 , is balanced by nonlinearity, represented by A. Here we 
have added to this balance weak inhomogeneity so that c x /c scales as h 2 d 3 /dx 3 , and weak 
friction so that Cd scales with hd/dx. Within this basic balance of terms, we can cast (JT]) 
into the asymptotically equivalent form 

^ + 'ik A + Yh AAx + ^™ = - c °w^ A ' < 2 > 

f x dx' 

where r = / , X = r — t. (3) 



jo c(0 

Here we have h = h(x(r)), explicitly dependent on the variable r which describes evolution 
along the path of the wave. 

The governing equation §2§ can be cast into several equivalent forms. That most com- 
monly used is the variable-coefficient KdV equation, obtained here by putting 

B = {ghf'^A (4) 

3 h g 1/A 

so that B T + 2gl/ih5/i BB x + ^ B xxx = ~C D jjj- A \B\B . (5) 

This form shows that, in the absence of friction term, i.e. when Cd = 0, equation (j2j) 
has two integrals of motion with the densities proportional to h l ^ A A and h 1 ! 2 A 2 . These 
are often referred to as laws for the conservation of "mass" and "momentum". However, 
these densities do not necessarily correspond to the corresponding physical entities. Indeed, 
to leading order, the "momentum" density is proportional to the wave action flux, while 
the "mass" density differs slightly from the actual mass density. This latter issue has been 
explored by Miles (1979), where it was shown that the difference is smaller than the error 
incurred in the derivation of equation and is due to reflected waves. 

Our main concern in this paper is with the behaviour of an undular bore over a slope 
in the presence of bottom friction, using the perturbed KdV equation (|2J), where we were 
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originally motivated by the possibility that the behaviour of a tsunami approaching the 
shore might be modeled in this way. The undular bore solution to the unperturbed KdV 
equation can be constructed using the well-known Gurevich-Pitaevskii (GP) (1974) approach 
(see also Fornberg and Whitham 1978). In this approach, the undular bore is represented 
as a modulated nonlinear periodic wave train. The main feature of this unsteady undular 
bore is the presence of a solitary wave (which is the limiting wave form of the periodic 
cnoidal wave) at its leading edge. The original initial-value problem for the KdV equation is 
then replaced by a certain boundary-value problem for the associated modulation Whitham 
equations. We note, however, that so far, the simplest, "(x/t)" -similarity solutions of the 
modulation equations have been used for the modelling of undular bores in various contexts 
(see Grimshaw and Smyth 1986, Smyth 1987 or Apel 2003 for instance). These solutions, 
while effectively describing many features of undular bores, are degenerate and fail to cap- 
ture, even qualitatively, some important effects associated with non-self-similar modulation 
dynamics. In particular, in the classical GP solution for the resolution of an initial jump 
in the unperturbed KdV equation, the amplitude of the lead solitary wave in the undular 
bore is constant (twice the value of the initial jump). On the other hand, the modulation 
solution for the undular bore evolving from a general monotonically decreasing initial profile 
shows that the lead solitary wave amplitude in fact grows with time (Gurevich, Krylov and 
Mazur 1989; Gurevich, Krylov and El 1992; Kamchatnov 2000). As we shall see, the very 
possibility of such variations in the modulated solutions of the unperturbed KdV equation 
has a very important fluid dynamics implication: in a general setting, the undular bore lead 
solitary wave cannot be treated as an individual KdV solitary wave but rather represents a 
part of the global nonlinear wave structure. In other words, while at every particular moment 
of time the lead solitary wave has the spatial profile of the familiar KdV soliton, generally, 
the temporal dependence of its amplitude cannot be obtained in the framework of single 
solitary wave perturbation theory. 

In the unperturbed KdV equation, the growth of the lead solitary wave amplitude is 
caused by the spatial inhomogeneity of the initial data. Here, however, the presence of a 
perturbation due to topography and/or friction serves as an alternative and/or additional 
cause for variation of the lead solitary wave amplitude. Thus, in the present case, the 
variation in the amplitude will have two components (which generally, of course, cannot be 
separated because of the nonlinear nature of the problem); one is local, described by the 
adiabatic perturbation theory for a single solitary wave, and the other one is nonlocal, which 
in principle requires the study of the full modulation solution. Depending on the relative 
values of the small parameters associated with the slope, friction and spatial non-uniformity 
of the initial modulations, we can take into account only one of these components, or a 
combination of them. 

The structure of the paper is as follows. First, in Section 2, we reformulate the basic model 
([T]) as a constant-coefficient KdV equation perturbed by terms representing topography 
and friction. Then we derive in Section 3 the associated perturbed Whitham modulation 
equations using methods recently developed by Kamchatnov (2004). Next, in Section 4, this 
Whitham system is integrated in the solitary- wave limit. Our purpose here is primarily to 
obtain the equation of a multiple characteristic, which defines the leading edge of a shoaling 
undular bore in the case when the modulations due to the combined action of the slope and 
bottom friction are small compared to the existing spatial modulations due to non-uniformity 
of the initial data. As a by-product of this integration, we reproduce and extend the known 
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results on the adiabatic variation of a single solitary wave (Miles 1983a,b). Then, in Section 
5, we carry out an analogous study of a cnoidal wave, propagating over a gradual slope and 
subject to friction, a case studied previously by Miles (1983 b) but under the restriction of 
zero mean flow, which is removed here. Finally, in Section 6 we study effects of a gradual 
slope and bottom friction on the front of an undular bore which represents a modulated 
cnoidal wave transforming into a system of weakly interacting solitons near its leading edge. 



2 Problem formulation 

For the purpose of the present paper it is convenient to recast (j2J) into the standard KdV 
equation form with constant coefficients, modified by certain perturbation terms. Thus we 
introduce the new variables 




so that Ut + QUUx + Uxxx = R = F(T)U - G(T)\U\U, (7) 



where F(T) = GiT) = AC D ^. (8) 

In this form, the governing equation (j7j) has the structure of the integrable KdV equation 
on the left-hand side, while the separate effects of the varying depth and the bottom friction 
are represented by the two terms on the right-hand side. This structure enables us to use 
the general theory developed in Kamchatnov (2004) for perturbed integrable systems. 

For much of the subsequent discussion, it is useful to assume that h(x) = constant, 
Cd = for x < in the original equation (pQ), which corresponds to F(T) = G(T) = for 
T < in (J7J). We shall also assume that A = for x > at t = 0, which corresponds to 
{7 = for X > on X = r{T) (see ([H])). Then we shall propose two types of initial- value 
problem for flTJ, and correspondingly for 

(a) Let a solitary wave of a given amplitude ao initially propagating over a flat bottom 
without friction (i.e a soliton described by an unperturbed KdV equation), enter the variable 
topography and bottom friction region at t = 0, x = (Fig. 1 a). 

(b) Let an undular bore of a given intensity propagate over a flat bottom without friction 
(the corresponding solution of the unperturbed KdV equation will be discussed in Section 
5). Let the lead solitary wave of this undular bore have the same amplitude ao and enter 
the variable topography and bottom friction region at t — 0, x — (Fig. lb). 

In particular, we shall be interested in the comparison of the slow evolution of these 
two, initially identical, solitary waves in the two different problems described above. The 
expected essential difference in the evolution is due to the fact that the lead solitary wave 
in the undular bore is generally not independent of the remaining part of the bore and can 
exhibit features that cannot be captured by a local perturbation analysis. The well-known 
example of such a behaviour, when a solitary wave is constrained by the condition of being 
a part of a global nonlinear wave structure, is provided by the undular bore solution of the 
KdV-Burgers (KdV-B) equation 

u t + Quu x + u xxx = fiu xx , /i < 1 . (9) 
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Figure 1: Isolated solitary wave (a) and undular bore (b) entering the variable topogra- 
phy/bottom friction region. 

Indeed, the undular bore solution of the KdV-B equation is known to have a solitary 
wave at its leading edge (see Johnson 1970; Gurevich & Pitaevskii 1987; Avilov, Krichever 
& Novikov 1987) and this solitary wave: (a) is asymptotically close to a soliton solution 
of the unperturbed KdV equation; and (b) has the amplitude, say a , that is constant in 
time. At the same time, it is clear that if one takes an isolated KdV soliton of the same 
amplitude a as initial data for the KdV-Burgers equation it would damp with time due 
to dissipation. The physical explanation of such a drastic difference in the behaviour of an 
isolated soliton and a lead solitary wave in the undular bore for the same weakly dissipative 
KdV-B equation is that the action of weak dissipation on an expanding undular bore is 
twofold: on the one hand, the dissipation tends to decrease the amplitude of the wave 
locally but, on the other hand, it "squeezes" the undular bore so that the interaction (i.e. 
momentum exchange) between separate solitons within the bore becomes stronger than in 
the absence of dissipation and this acts as the amplitude increasing factor. The additional 
momentum is extracted from the upstream flow with a greater depth (see Benjamin and 
Lighthill 1954). As a result, in the case of the KdV-B equation, an equilibrium non-zero 
value for the lead solitary wave amplitude in the undular bore is established. Of course, 
for other types of dissipation, a stationary value of the lead soliton amplitude would not 
necessarily exist, but in general, due to the expected increase of the soliton interactions near 
the leading edge, the amplitude of the lead soliton of the undular bore would decay slower 
than that of an isolated soliton. Indeed, the presence here of variable topography as well 
can result in an additional "nonlocal" amplitude growth. 

While the problem (a) can be solved using traditional perturbation analysis for a single 
solitary wave, which leads to an ordinary differential equation along the solitary wave path 
(see Miles 1983a,b), the undular bore evolution problem (b) requires a more general approach 
which can be developed on the basis of Whitham's modulation theory leading to a system of 
three nonlinear hyperbolic partial differential equations of the first order. Since the Whitham 



5 



method, being equivalent to a nonlinear multiple scale perturbation procedure, contains the 
adiabatic theory of slow evolution of a single solitary wave as a particular (albeit singular) 
limit, it is instructive for the purposes of this paper to treat both problems (a) and (b) using 
the general Whitham theory. 



3 Modulation equations 

The original Whitham method (Whitham 1965, 1974) was developed for conservative constant- 
coefficient nonlinear dispersive equations and is based on the averaging of appropriate con- 
servation laws of the original system over the period of a single-phase periodic travelling 
wave solution. The resulting system of quasi-linear equations describes the slow evolution 
of the modulations (i.e. of the mean value, the wavenumber, the amplitude etc.) of the pe- 
riodic travelling wave. Here, that approach is extended to the perturbed KdV equation (jBJ) 
following the general approach of Kamchatnov (2004), which extends earlier results for cer- 
tain specific cases (see Gurevich and Pitaevskii (1987, 1991), Avilov, Krichever and Novikov 
(1987) and Myint and Grimshaw (1995) for instance). 

We suppose that the evolution of the nonlinear wave is adiabatically slow, that is, the 
wave can be locally represented as a solution of the corresponding unperturbed KdV equation 
(i.e. (j7]) with zero on the right-hand side) with its parameters slowly varying with space and 
time. The one-phase periodic solution of the KdV equation can be written in the form 



U(X,T) = A 3 - Ax - A 2 - 2(A 3 - A 2 )sn 2 (v / A 3 - Xi8,m) (10) 

where sn(y, m) is the Jacobi elliptic sine function, Ai < A 2 < A3 are parameters and the 
phase variable 9 and the modulus m are given by 

9 = X-VT, V = -2(A! + A 2 + A 3 ), (11) 
^3 - A 2 

m=- -, (12) 



A3 — Ai 



a 3 



dfi 2K(m) 



and L= <bd9= . = = . v ' , (13) 



A 2 



where K(m) is the complete elliptic integral of the first kind, L is the "wavelength" along the 
X-axis (which is actually a retarded time rather than a true spatial co-ordinate). Here we 
have used the representation of the basic ordinary differential equation for the KdV travelling 
wave solution ({TO]) in the form (see Kamchatnov (2000) for a general motivation behind this 
representation) 



where 
and 



H = \(U + s x ), si = A! + A 2 + A 3 (15) 

3 
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that is the solution ([TO]) is parameterized by the zeroes Ai, A 2 , A 3 of the polynomial 

In a modulated wave, the parameters Ai, A 2 , A3 are allowed to be slow functions of X and 
T, and their evolution is governed by the Whitham equations. For the unperturbed KdV 
equation, the evolution of the modulation parameters is due to a spatial non-uniformity 
of the initial distributions for Xj,j = 1,2,3 and the typical spatio-temporal scale of the 
modulation variations is determined by the scale of the initial data. 

In the case of the perturbed KdV equation (CO), the evolution of the parameters Ai, A 2 , A3 
is caused not only by their initial spatial non-uniformity, but also by the action of the 
weak perturbation, so that, generally, at least two independent spatio-temporal scales for 
the modulations can be involved. However, at this point we shall not introduce any scale 
separation within the modulation theory and derive general perturbed Whitham equations 
assuming that the typical values of F(T) and G(T) are 0(dXj/dT,dXj/dX) within the 
modulation theory. 

It is instructive to first introduce the Whitham equations for the perturbed KdV equation 
(J?D using the traditional approach of averaging the (perturbed) conservation laws. To this 
end, we introduce the averaging over the period fflBl of the cnoidal wave ffTU]) by 



In particular, 



(U) = 2(ji) - 8l = 2(A 3 - Xi)^r\ + Ai - A 2 - A 3 , 

K{m) 

T2\ _ or s l/\ x ^ E ( m ) * \ , !/\2 



(U 2 ) = 8[-^(A 3 - Ai)-^ - - ai Ai + -(A - A 2 A 3 )] + s , (19) 
K{m) 3 

where E(m) is the complete elliptic integral of the second kind. Now, one represents the 
KdV equation ([7j) in the form of the perturbed conservation laws 



§T + 1I = jRj ' j L ' 2 ' 3 ' R > < : m 



where Pj and Qj are the standard expressions for the conserved densities (Kruskal integrals) 
and "fluxes" of the unperturbed KdV equation. Just as in the Whitham (1965) theory for 
unperturbed dispersive systems, the number of conservation laws required is equal to the 
number of free parameters in the travelling wave solution, which is three in the present case. 
Next, one applies the averaging fTTTj) to the system fT20l) to obtain (see Dubrovin and Novikov 
1989) 

d(P s ) d{Q j )_ 

~df~ + ~dX~ ~ ^ R i> ' •7- 1 > 2 ' 3 - ( 21 ) 

The system (12ip describes slow evolution of the parameters Xj in the cnoidal wave solution 

Along with these derived perturbed conservative form of the Whitham equations, we 
introduce the wave conservation law which is a general condition for the existence of slowly 
modulated single-phase travelling wave solutions ([TO]) (see for instance Whitham 1974) and 
must be consistent with the modulation system ( l2~li) . This conservation law has the form 

dk duo , . 

9T + ax=°- < 22 > 
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where k = — , uj = kV (23) 
lj 

are the "wavenumber" and the "frequency" respectively (we have put quotation marks here 
because the actual wavenumber and frequency related to the physical variables x, t are 
different quantities from those in ( l23l) . but are related through the transformations (jSlEl) )• 
The wave conservation law ( l22i) can be introduced instead of any of three inhomogeneous 
averaged conservation laws comprising the Whitham system ( T2T1 ). 

It is known that the Whitham system for the homogeneous constant-coefficient KdV 
equation can be represented in diagonal (Riemann) form (Whitham 1965, 1974) by an ap- 
propriate choice of the three parameters characterising the periodic travelling wave solution. 
In fact, in our solution ([7]) the parameters Xj have already been chosen so that they coincide 
with the Riemann invariants of the unperturbed KdV modulation system. Introducing them 
explicitly into the perturbed system (1211) we obtain (see Kamchatnov 2004) 

8T +V% dX- dL/dXi 4n^(A,-A,) ' 1 ~ ^ 
where R is the perturbation term on the right-hand side of the KdV equation ([7j) and 

are the Whitham characteristic velocities corresponding to the unperturbed KdV equation. 

It should be noted that the straightforward realisation of the above lucid general algo- 
rithm for obtaining perturbed modulation system in diagonal form is quite a laborious task. 
In fact, to derive system ( 124"1) . the so-called finite-gap integration method incorporating the 
integrable structure of the unperturbed KdV equation has been used. The modulation sys- 
tem (1241) in a more particular form corresponding to specific choices of the perturbation term 
was obtained by Myint and Grimshaw (1995) using a multiple-scale perturbation expansion. 
In that latter setting, the wave conservation law f[2"2"j) is an inherent part of the construction, 
while in the averaging approach used here, it can be obtained as a consequence of the system 

dUD- 

To obtain an explicit representation of the Whitham equations for the present case of 
equation we must substitute the perturbation R from the right-hand side of © and 
perform the integration 0171) with U given by (11 01) . From now on, we are going to consider 
only the flows where U > so that the perturbation term assumes the form 

R(U) =G(T)U - F(T)U 2 . (26) 

Substituting (|2"6"|) into §M§ we obtain, after some detailed calculations (see Appendix), 
the perturbed Whitham system in the form 

^ + "'^ = ft = Ci[F(m ~ G(r,Bi1, s = 1 ' 2 - 3 < 27 > 

where Ci = i c >= E-(i- m) K > c > = Eh<-- (28 > 



s 



At = i(5Ai - A 2 - X 3 )E + jj(A 2 - Xt)K, 

A 2 = ^(5A 2 - Ai - A 3 )£ - (A 2 - Ai) Q + ) /f, ( 29 ) 

if; 



A 3 = i(5A3-A!-A 2 )E- 



A 3 + -(A 2 -Ai) 



Si = ^r(-27A? - 7A^ - 7A" + 2A X A 2 + 2AiA 3 + 22A 2 A 3 )£ 
15 

-^-(\2-Xi)(SX 1 + X 2 + X 3 )K, 
15 

B 2 = ^r(-7A? - 27A" - 7Xj + 2AxA 2 + 22AiA 3 + 2X 2 X 3 )E 

+ 4 a 2 " A 1 ^! + 15A 2 + 11A 3 - 6A 1 A 2 - 18A 1 A 3 + 6A 2 A 3 )iT, 

15 A 3 — Ai 

B 3 = j^(-7Xl - 1X\ - 11X1 + 22A i A 2 + 2AxA 3 + 2A 2 A 3 )£ 

+ ^r(7A? + 11A^ + 15A* - 18AiA 2 - 6AiA 3 + 6A 2 A 3 )if; 
and the characteristic velocities are: 

4(A 3 -Ai)(l-m)tf 



(30) 



^3 = -2 ^ A i 



4(A 3 - A 2 )AT 
E — K ' 



The equations (12"T|) - (13"Tj) provide a general setting for studying the nonlinear modulated 
wave evolution over variable topography with bottom friction. In the absence of the pertur- 
bation terms (i.e. when F(T) = 0, G{T) = 0), the system (1271) . (1311) indeed coincides with 
the original Whitham equations (Whitham 1965) for the integrable KdV dynamics. In that 
case the variables Aj, A 2 , A 3 become Riemann invariants, so in this general (perturbed) case 
we shall call them Riemann variables. 

It is important to study the structure of the perturbed Whitham equations (1271) - (I3TI) in 
two limiting cases when the underlying cnoidal wave degenerates into (i) a small-amplitude 
sinusoidal wave (linear limit), when A 2 = A 3 (m = 0), and (ii) into a solitary wave when 
X 2 = Ai (m = 1). Since in both these limits the oscillations do not contribute to the mean 
flow (they are infinitely small in the linear limit and the distance between them becomes 
infinitely long in the solitary wave limit) one should expect that in both cases one of the 
Whitham equations will transform into the dispersionless limit of the original perturbed 
KdV equation (jTJ) i.e. 

U T + 6UU X = F{T)U - G(T)U 2 , (32) 
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Indeed, using formulae (127|) - (!3T|) we obtain for m = 0: 



A2 — A3 , 

~ 6Al 9X =AlF + AlG ' ( 33 ) 
^ + (6A 1 -12A 3 )^ = A 1 F + A?G. 

Similarly, for m = 1, one has 
A2 = Ai , 

<9A <9A 1 1 

_1 _ ( 4Al + 2A 3 )^ = -(4A, - A 3 )F + -(7A| - 24A t A 3 + 32A?)G, (34) 

We see that, in both cases, one of the Riemann variables (taken with inverted sign) coincides 
with the solution of the dispersionless equation fl32l) (recall that in the derivation of the 
Whitham equations we assumed U > everywhere), namely U = (U) = — Ai when A 2 = A 3 
(m = 0) and U = (U) = —A3 when A2 = Ai (m = 1). 

To conclude this section, we present expressions for the physical wave parameters such as 
the surface elevation wave amplitude a, mean elevation (A) speed and wavenumber in terms 
of the modulation solution Xj(X, T). Using (J6]) and (jTUI) we obtain for the wave amplitude 
(peak to trough) and the mean elevation 

Ah 2 2h 2 
a = —(X 3 -X 2 ), ( A ) = ^(U), (35) 

where the dependence of (U) on Xj(X,T), j = 1,2,3 is given by f|T8|) and X = X(x,t), 
T = T(x,t) by ([31 [6]). In order to obtain the physical wavenumber k and the frequency Q 
we first note that the phase function 6(X,T) defined in ffTTl) is replaced by a more general 
expression defined so that k = 9x and kV = —6? are the "wavenumber" and "frequency" in 
the X — T coordinate system. Then we define the physical phase function 0(x, t) = 6(X, T) 
so that we get 

k = b x , n = -e t . (36) 



It now follows that 



k^ i hV , „ _ VL c , _. 

, Q = k, and - = - — — . (37) 

c bg k 1 — liv/bg 



Note that the physical frequency is the "wavenumber" in the X — T coordinate system, 
and that the physical phase speed is VL/k. Since the validity of the KdV model (p]) requires 
inter alia that the wave be right-going it follows from this expression that the modulation 
solution remains valid only when hV < bg. Of course, the validity of (CQ) also requires that 
the amplitude remains small, and this would normally also ensure that V remains small. 
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4 Modulation solution in the solitary wave limit 



In this section, we shall integrate the perturbed modulation system fl27j) along the multiple 
characteristic corresponding to the merging of two Riemann variables A2 and Ai. As we shall 
see later, this characteristic specifies the motion of the leading edge of the shoaling undular 
bore in the case when the perturbations due to variable topography and bottom friction can 
be considered as small compared with the existing spatial modulations within the bore. At 
the same time, as the case A2 = Ai ( i.e. m — 1) corresponds to the solitary wave limit in the 
travelling wave solution fflQj) . our results here are expected to be consistent with the results 
from the traditional perturbation approach to the adiabatic variation of a solitary wave due 
to topography and bottom friction (see Miles 1983a,b). 

In the limit m — > 1 the periodic solution ( flOj) of the KdV equation goes over to its solitary 
wave solution 



are the solitary wave amplitude and "velocity" respectively. The solution (1381) depends on two 
parameters Ai and A3 whose adiabatic slow evolution is governed by the reduced modulation 
system fl34|) . It is important that the second equation in this system is decoupled from 
the first one. Hence, evolution of the pedestal —A3 on which the solitary wave rides, can 
be found from the solution of this dispersionless equation by the method of characteristics. 
When A3(A, T) is known, evolution of the parameter Ai can be found from the solution of 
the first equation f[3~4l) . As a result, we arrive at a complete description of adiabatic slow 
evolution of the solitary wave parameters taking account of its interaction with the (given) 
pedestal. 

However, it is important to note here that while this description of the adiabatic evolution 
of a solitary wave is complete as far as the solitary wave itself is concerned, it fails to describe 
the evolution of a trailing shelf, which is needed to conserve total "mass" (see, for instance, 
Johnson 1973b, Grimshaw 1979 or Grimshaw 2006). This trailing shelf has a very small 
amplitude, but a very large length scale, and hence can carry the same order of "mass" as 
the solitary wave. But note that the "momentum" of the trailing shelf is much smaller than 
that of the solitary wave, whose adiabatic deformation is in fact governed to leading order by 
conservation of "momentum" , or more precisely, by conservation of wave action flux (strictly 
speaking, conservation only in the absence of friction). 

The situation simplifies if the solitary wave propagates into a region of still water so 
that there is no pedestal ahead of the wave, that is A3 = in X > r(T). But then, since 
A3 = is an exact solution of the degenerate Whitham system (1341) for this solitary wave 
configuration, we can put A3 = both in the solitary wave solution, 




(38) 



where 




(39) 



U{X,T) = -2Aisech 2 [^A^pf - V S T)), V s = -4A 



(40) 



and in equation ( 134"1) for the parameter Ai to obtain, 



1 



-4A 



dX 




(41) 



dT 



OX 



1 ■ 
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As we see, the solitary wave moves with the instant velocity 

dX 
~dT 



-4Ai, 



(42) 



and the parameter Ai changes with T along the solitary wave trajectory according to the 
ordinary differential equation 



d\-\ 4 , 32 , n 



(43) 



It can be shown that equation ( 1431) is consistent with the equation for the solitary wave half- 
width 7 = \J — \\ obtained by the traditional perturbation approach (see Grimshaw (1979) 
for instance). 

Next, we re- write equation ([43]) in terms the original independent x- variable. For that, 
we find from ([6]), that 

dT = {h^/Qg^dx (44) 

21 f9 ^dh G = AC f9\^ 
dx \h, 



and F 



- 



2 \hJ 

Then substituting these expressions into (j4"3j) yields the equation 

d\i ^ 1 dh 
h dx 



(45) 



dx 



_ 3f = Al + 6 4^M 



which can be easily integrated to give 



45 g 



64 C D [ x dx 



(46) 



(47) 



where Co is an integration constant and x = is a reference point where h = ho- According 
to ( )40|) . £/"o = — 2Ai is the amplitude of the soliton expressed in terms of variable U(X, T). 
Returning to the original surface displacement A(x, t) by means of and denoting Co = 
4/(3gaoh ), we find the dependence of the surface elevation soliton amplitude a = (2h 2 /3g)Uo 
on x in the form 



a 



ho 
h 



1 + —C D a ho 
15 



x dx 



-i 



(48) 



where a>o is the solitary wave amplitude at x — 0. We note that for Cp = this reduces to 
the classical Boussinesq (1872) result a ~ h" 1 , while for h = ho it reduces to the well-known 
algebraic decay law a ~ 1/(1 + constant x) due to Chezy friction. Miles (1983a,b) obtained 
this expression for a linear depth variation, although we note that there is a factor of 2 
difference from (HHj) (in Miles (1983a,b) the factor lQCn/15 is 8Cx>/15). The trajectory of 
the soliton can be now found from (H2|) and fj4"Tj) : 



A 



aoho 



dx'h- 5 / 2 (x') 



16 
15 



Cndoho 



dx 
h s (x) 



(49) 



This expression determines implicitly the dependence of x on t along the solitary wave path 
and provides the desired equation for the multiple characteristic of the modulation system 
for the case m — 1. 
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It is instructive to derive an explicit expression for the solitary wave speed by computing 
the derivative dx/dt from (1491) . or more simply, directly from (1371) . 



dx 
~dt 



a/2h 



(50) 



The formula fl50l) yields the restriction for the relative amplitude 7 = a/h < 2 which is 
clearly beyond the applicability of the KdV approximation (wave breaking occurs already 
at 7 = 0.7 (see Whitham 1974)). In the frictionless case {Cp = 0) equation (T4"51) gives 
a/h = aoh Q /h 2 , and so the expression (15"U1) for the speed must fail as h —>■ 0. It is interesting 
to note that this failure of the KdV model as h — > due to appearance of infinite (and 
further negative!) solitary wave speeds is not apparent from the expression (148]) for the 
solitary wave amplitude, and the implication is that the model cannot be continued as 
h — > 0. Curiously this restriction of the KdV model seems never to have been noticed before 
in spite of numerous works on this subject. Note that taking account of bottom friction 
leads to a more complicated formula for the solitary wave speed as a function of h but the 
qualitative result remains the same. 

It is straightforward to show from fj46l) or fT48l) that 



Ox 

a 



h x lQC D a o h 



h 15 h 3 



, j [ x dx 

1 + T^^Daoho / — 



15 



h 3 



(51) 



It follows immediately that for a wave advancing into increasing depth (h x > 0), the ampli- 
tude decreases due to a combination of increasing depth and bottom friction. However, for 
a wave advancing into decreasing depth, there is a tendency to increase the amplitude due 
to the depth decrease, but to decrease the amplitude due to bottom friction. Hence whether 
or not the amplitude increases is determined by which of these effects is larger, and this in 
turn is determined by the slope, the depth, and the consolidated drag parameter Coao/ho- 
To illustrate, let us consider the bottom topography in the form 



h(x) 



hl- a (h -5xY 



a > 



(52) 



which satisfies the condition h(0) = ho] the parameter 5 characterizes the slope of the bot- 
tom. In this case the formula (1481) becomes 



h 



1 + 



16 Codo 
15 5(3a- l)ho 




1 



(53) 



if a 7^ 1/3. One can see now that if a < 1/3, then the bottom friction term is relatively 
unimportant due to the smallness of Cu- Of course, for this case we again recover the 
Boussinesq result, now slightly modified, 



a ~ a - 



ho 
h 



1 + 



16 Cdclo 
15 5(1 -3a)h. 



n -1 



< a < 



1 



h < h c 



(54) 



Of course, this result is impractical in the KdV context as the KdV approximation used here 
requires the ratio a/h to remain small. 
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If a > 1/3 now obtain asymptotic formula 



15(3a- 1)5 



i-2 



/l< /in 



(55) 



which is independent of the initial amplitude an. This expression is consistent with the small- 
amplitude KdV approximation as long as (3a — 1)5 /Co is order unity. Simple inspection of 
( |55i) shows that the solitary wave amplitude 



increases as h 



if | < a < §, 



• is constant as h —>■ if a = |, 

• decreases as/i— ► if a > |. 

Thus for 1/3 < a < 1/2, as for the case a < 1/3, the amplitude will increase as the depth 
decreases, in spite of the presence of (sufficiently small) friction. However, for a > 1/3, even 
although there is usually some initial growth in the amplitude, eventually even small bottom 
friction will take effect and the amplitude decreases to zero. We note that if a = 1/3 then 
the integral f Q h~ 3 dx in (j4"8]) diverges logarithmically as h — > 0, which just slightly modifies 
the result (1551) for h <C h Q and implies growth of the amplitude oc In h/h as h — > 0. 
Of particular interest is the case a = 1. In that case formula (1531 becomes 



a 



ho 
h 

and 




-i 



15 5 , 



h < h 



(56) 



(57) 



These expressions (|56l [57|) were obtained by Miles (1983a,b) using wave energy conservation 
(as above, note, however, that in Miles (1983a,b) the numerical coefficient is 15/4 rather 
than 15/8). Thus, these results obtained from the Whitham theory are indeed consistent, at 
the leading order, with the traditional perturbation approach for a slowly-varying solitary 
wave. 



5 Adiabatic deformation of a cnoidal wave 

Next we consider a modulated cnoidal wave fflOl) in the special case when the modulation does 
not depend on X. While this case is, strictly speaking, impractical as it assumes there is an 
infinitely long wavetrain, it can nevertheless provide some useful insights into the qualitative 
effects of gradual slope and friction on undular bores which are locally represented as cnoidal 
waves. In the absence of friction, the slow dependence of the cnoidal wave parameters on T 
was obtained by Ostrovsky & Pelinovsky (1970, 1975) and Miles (1979) (see also Grimshaw 
2006), assuming that the surface displacement had a zero mean (i.e. (U) = 0), while, 
the effects of friction were taken into account by Miles (1983b) using the same zero- mean 
displacement assumption. However, this assumption is inconsistent with our aim to study 
undular bores where the value of (U) is essentially nonzero. Hence, we need to develop a 
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more general theory enabling us to take into account variations in all the parameters in the 
cnoidal wave. Such a general setting is provided by the modulation system (j27j) . 

Thus we consider the case when the Riemann variables in (1271) do not depend on the 
variable X so that the general Whitham equations become ordinary differential equations 
in T, which can be conveniently reformulated in terms of the original spatial x-coordinate 
using the relationship 



— =C- 

dx 



9 1dh A 2C D 
4hdx 1 3g % 



1,2,3, (5* 



where all variables are defined above in section 3 (seeEElEHlEO]). This system can be readily 
solved numerically. But it is instructive, however, to first indicate some general properties 
of the solution. 

First, the solution to the system (!58|) must have the property of conservation of "wave- 
length" L (or "wavenumber" k=27r/L) 

2K(m) . . 

L = = = constant (59) 

V A3 — Ai 

Indeed, the wave conservation law ( 1221) in absence of X-dependence assumes the form 

dk , . 

which yields ( 159]) . Thus the system of three equations (!58j) can be reduced to two equations. 
Next, applying Whitham averaging directly to ([7]) yields 

™ -HfM-?^P, M = (U) , P= m u). (61) 
dx Andx Sg 

The equation set (159~j) . (ISTj) . (1521 comprise a closed modulation system for three independent 
modulation parameters, say M, P and m. While this system is not as convenient for further 
analysis as the system (1271 in Riemann variables, it does not have a restriction U > inherent 
in ( 127]) . and allows for some straightforward inferences regarding the possible existence of 
modulation solutions with zero mean elevation, that is with M — 0. Indeed, one can see 
that the solution with the zero mean is actually not generally permissible when Cp 7^ 0, a 
situation overlooked in Miles (1983b). Indeed, M = immediately then implies that P = 
by (16T1) . But then due to (1591) we have all three modulation parameters fixed which is clearly 
inconsistent with the remaining equation (162|) (except for the trivial case M = 0, P = 0, 
Q — 0). However, in the absence of friction, when Co = 0, equation (16T|) uncouples and 
permits a nontrivial solution with a zero mean. In general, when Co = equations (1610 . 
(16"2|) can be easily integrated to give 

d = Mh 9/4 = constant; cr = Ph 9/2 = constant. (63) 

Then, using ( 118] [T91 15^]) one readily gets the formula for the variation of the modulus m, 
and hence of all the other wave parameters, as a function of h 

K\2{2 - m)EK - 3£ 2 - (1 - m)K>] = " { ° "ff L * . (64) 
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m 




Figure 2: Dependence of the modulus m on the physical space coordinate x in the cases 
without and with bottom friction in the AT-mdependent modulation solution. 



Formula (164p generalises to the case M ^ (i.e. d ^ 0) the expressions of Ostrovsky & 
Pelinovsky (1970, 1975), Miles (1979) and Grimshaw (2006) (note that in Grimshaw (2006) 
the zero mean restriction in actually not necessary). We note here that, again with Co = 0, 
equation ([5]) implies conservation of (B) and (B 2 ) (the averaged wave action flux), which, 
together with (|59l . also yield (|64l) . 

The physical frequency Q and wavenumber k in the modulated periodic wave under study 
are given by the formula ( 1371) . and we recall here that k = 2n / L is constant (see (1591) ). As 
discussed before at the end of Section 3 we must require that the phase speed stays positive as 
the wave evolves, and here that requires that the physical wavenumber k > 0. Since a/h (and 
hence hV/6g) is supposed to be small within the range of applicability of the KdV equation 
(T5]) the expression ( 1371) implies the behaviour k ~ Q/y/gh which of course agrees with the 
well known result for linear waves on a sloping beach (see Johnson 1997 for instance). This 
effect will be slightly attenuated for the nonlinear cnoidal wave, since Vh/6g > 0, but the 
overall effect will be a "squeezing" of the cnoidal wave, a result important for our further 
study of undular bores. Next we study numerically the combined effect of slope and friction 
on a cnoidal wave. 

As we have shown, in the presence of Chezy friction M^fl, and we have also assumed 
that U > 0, which is necessary when we come to study undular bores. Now we use the 
stationary modulation system (158|) in Riemann variables, which was derived using this as- 
sumption. We solve the coupled ordinary differential equation system (158|) for the case of a 
linear slope 

h(x) = h — 5x (65) 
with ho = 10, 5 = 0.01, and with the initial conditions 

Ai = -0.441, A 2 = 0.147, A 3 = 0.294 at x = 0, (66) 

which corresponds to a nearly harmonic wave with m = 0.2, a/ ho = 0.2, (A) /ho w 0.3 
at x = (see (135]) ). Also we note that for the chosen parameters we have V — 0, so at 
x = we have k = Q/^gho as in linear theory It is instructive to compare solutions with 
(Co = 0.01) and without (Co = 0) friction. In Fig. [2] the dependence of the modulus m 
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<A> 



h 1/4 <A> 



c D= o 




Figure 3: Left: Dependence of the mean value (A) in the X-independent modulation solution 
on the physical space coordinate x without (dashed line) and with (solid line) bottom friction; 
Right: Same but multiplied by the Green's law factor, h 1 ^ 



C D =0 




Figure 4: Dependence of the surface elevation amplitude a on the space coordinate x. Dashed 
line corresponds to the frictionless case and solid line to the case with bottom friction. 



on x is shown for both cases. We see that for the frictionless case m — > 1 with decrease 
of depth, i.e. the wave crests assume the shape of solitary waves when one approaches the 
shoreline. When Co ^ the modulus also grows with decrease of depth but never reaches 
unity. The dependence on x of the mean surface elevation (A) for the cases without and 
with friction is shown in Fig. [3j We have checked that the "wavelength" L (1591) is constant 
for both solutions. Also, one can see from Fig. [3] (right) that the value /i 1 / 4 (A) oc d is indeed 
conserved in the frictionless case but is not constant if friction is present (the same holds 
true for the value ^^(A 2 ) oc a but we do not present the graph here). Finally, in Fig. H] the 
dependence of the physical elevation wave amplitude a on the spatial coordinate x is shown. 
One can see that the amplitude adiabatically grows with distance in the frictionless case 
due to the effect of the slope (without friction) but, not unexpectedly, gradually decreases 
in the case when bottom friction is present, where the decrease for these parameter settings 
is comparable in magnitude to the effect of the slope. In both cases the main qualitative 
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changes occur in the wave shape and the wavelength. 

Overall, we can infer from these results that the main local effect of a slope and bottom 
friction on a cnoidal wave, along with the adiabatic amplitude variations, is twofold: a wave 
with a m < 1 at j; = tends to transform into a sequence of solitary waves as x decreases, 
and at the same time the distance between subsequent wave crests tends to decrease. This 
is in sharp contrast with the behaviour of modulated cnoidal waves in problems described 
by the unperturbed KdV equation, where growth of the modulus m is accompanied by an 
increase of the distance between the wave crests. Generally, in the study of behaviour of 
unsteady undular bores in the presence of a slope and bottom friction we will have to deal 
with the combination of these two opposite tendencies. 

6 Undular bore propagation over variable topography 
with bottom friction 



6.1 Gurevich-Pitaevskii problem for flat-bottom zero-friction case 

We now turn to the problem (b) outlined in Section 2. We study the evolution of an undular 
bore developing from an initial surface elevation jump A > 0, located at some point xq < 0. 
As discussed below, the undular bore will expand with time so that at some t — to its lead 
solitary wave enters the gradual slope region, which begins at x = (see Fig. lb). We assume 
that for x < one has h = h Q = constant and Cd = 0. We shall first present a formulation 
of the Gurevich-Pitaevskii problem for the perturbation-free KdV equation and reproduce 
the well-known similarity modulation solution describing the evolution of the undular bore 
until the moment it enters the slope. We emphasize that, although this formulation and, 
especially, this similarity solution are known very well and have been used by many authors, 
some of the inferences important for the present application to fluid dynamics have not been 
widely appreciated, as far as we can discern. Pertinent to our main objective in this paper, 
we undertake a detailed study of the characteristics of the Whitham modulation system in 
the vicinity of the leading edge of the undular bore solution, and show that the boundary con- 
ditions of Gurevich-Pitaevskii type permit only two possible characteristics configurations, 
implying two qualitatively different types of the leading solitary wave behaviour. Next, we 
shall show how this Gurevich-Pitaevskii formulation of the problem applies to the perturbed 
modulation system in the form fTSTl) and finally we will study the effects of the perturbation 
on the modulations in the vicinity of the leading edge of the undular bore. 

In the case of a flat, frictionless bottom the original equation ([I]) becomes the constant- 
coefficient KdV equation which can be cast into the standard form 

77c + Q VVz + = (67) 

by introducing the new variables 

where Xq < is an arbitrary constant. In the Gurevich-Pitaevskii (GP) approach, one 
considers a large-scale initial disturbance ?7(£,0) = /(C), in the form of a decreasing profile, 
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/'(£) < (e.g. a smooth step: /(£) — >■ as £ — > +00; /(£) — > r] > as £ — > —00), whose 
initial evolution until some critical (breaking) time (b can be described by the dispersionless 
limit of the KdV equation, i.e. by the Hopf equation, 

(<( b : »j«r(£,0, r c + 6rr ? = 0, r(£, 0) = /(£) . (69) 

The evolution ( 1691) leads to wave-breaking of the r(£)-profile at some £ = with the 
consequence that the dispersive term in the KdV equation then comes into play, and an 
undular bore forms, which can be locally represented as a single-phase travelling wave. This 
travelling wave is modulated in such a way that it acquires the form of a solitary wave at the 
leading edge £ = £ + (C) an d gradually degenerates, via the nonlinear cnoidal-wave regime, to 
a linear wave packet at the trailing edge £ = £~(C)- It is important that this undular bore 
is essentially unsteady, i.e. the region £~(£) < £ < £ + (C) expands with time £. 

The single-phase travelling wave solution of the KdV equation (l6"7j) has the form (cf. 

(tH) 

V{€> = r 3 - n - r 2 - 2(r 3 - r 2 )sn 2 ( v / r 3 - n9, m) (70) 

fl = £ + 2(r 1 + r 2 + r 3 )C, m = . (71) 

r 3 - n 

The parameters rj < r 2 < r 3 < in the undular bore are slowly varying functions of £, (, 
whose evolution is governed by the Whitham equations 

OT ■ OT ■ 

-jr + v j (r 1 ,r 2 ,r 3 )-± = 0, 3 = 1,2,3. (72) 

The characteristic velocities in (1721 are given by (I3TI) . We stress that, although analytical 
expressions (170]) and (ITU|) (as well as (1721 and the homogeneous version of (12"7I) ) are identical, 
they are written for completely different sets of variables, both dependent and independent. 

The Riemann invariants rj(£, £) are subject to special matching conditions at the free 
boundaries, £ = £ ± (C) defined by the conditions m = (trailing edge) and m = 1 (leading 
edge), formulated in Gurevich and Pitaevskii (1974) (see also Kamchatnov (2000) or El 
(2005) for a detailed description). 

At the trailing (harmonic) edge, where the wave amplitude a = 2(r 3 — r 2 ) vanishes and 
m = 0, one has 

£ = <T(C): r 2 = r 3 , -r x = r. (73) 
At the leading (soliton) edge, where m = 1 one has 

£ = £ + (C): r 2 = r 1; -r 3 = r. (74) 



In both (1731 and (1741 . r(£, £) is the solution of the Hopf equation 

The curves £ = £±(() are defined for the solution of the GP problem (172|) . (173|) . (174j) by 
the ordinary differential equations 

£ = .-(r,o, f- + « + .o, (75) 

where t> ± are calculated as the values of double characteristic velocities of the modulation 
system at the undular bore edges, 

V = v 2 (r h r 3 ,r 3 )\^- {0 = v 3 (ri,r 3 ,r 3 )\^- {0 , (76) 



19 



v + = v 2 {r 1 ,r 1 ,r 3 )\^+ {0 = v x (r u r u r 3 )| €=c+(c) (77) 

These equations ( 1751) essentially represent kinematic boundary conditions for the undular 
bore (see El 2005). Indeed, the double characteristic velocity V2(ri,r 3 ,r 3 ) = v 3 (ri, r 3 ,r 3 ) can 
be shown to coincide with the linear group velocity of the small-amplitude KdV wavepacket 
while the double characteristic velocity ^2(^1,^1,^3) = Vi(ri, r±, r 3 ) is the soliton speed. 

One might infer from this GP formulation of the problem that, since the leading edge of 
the undular bore specified by (175]) . (1771) is a characteristic of the modulation system, then 
the value of the double Riemann invariant r + = r 2 = r\ is constant. Then, on considering an 
undular bore propagating into still water, where r = 0, one would obtain from the matching 
condition ([71]) at the leading edge that r 3 |^ = ^+ = and thus, the amplitude of the lead solitary 
wave a + = 2(r 3 — r 1 )\^+ = —r + would always be constant as well. However, this contradicts 
the general physical reasoning that the amplitude of the lead solitary wave should be allowed 
to change in the case of general initial data. The apparent contradiction is resolved by noting 
that the leading edge specified by (1751) , ( 1771) can be an envelope of the characteristic family, 
i.e. a caustic, rather than necessarily a regular characteristic, and hence there is no necessity 
for the double Riemann invariant r + to be constant along the curve £ = £ + (C) m general case. 
On the other hand, since the leading edge is defined by the condition m — 1, the wave form 
at the leading edge will coincide with the spatial profile of the standard KdV soliton. Thus 
we arrive at the conclusion that, in general, the amplitude of the leading KdV solitary wave 
will vary, even in the absence of the perturbation terms. Of course, in the unperturbed KdV 
equation, such varying solitary waves cannot not exist on their own, and require the presence 
of the rest of the undular bore. We also stress that these variations of the leading solitary 
wave in the undular bore, as described here, have a completely different physical nature to 
the variations of the parameters of an individual solitary wave due to small perturbations as 
described in Section 4. They are caused by nonlinear wave interactions within the undular 
bore rather than by a local adiabatic response of the solitary wave to a perturbation induced 
by topography and friction. Importantly for our study, however, it will transpire that the 
action of these same perturbation terms on the undular bore can lead to both a local and a 
nonlocal response of the leading solitary wave. 

6.2 Undular bore developing from an initial jump 

Next we consider the simplest solution of the modulation system, which describes an undular 
bore developing from an initial discontinuity placed at the point x = —Xq. In (77; £, £) - 
variables we have the initial conditions 

7?(£,0) = A for £<0; 77(C, 0) = for £>0, (78) 

where A > is a constant. Then, on using (|69|) . the initial conditions (178]) are readily 
translated into the free-boundary matching conditions (!73|) . (1741) for the Riemann invariants. 
Because of the absence of a length scale in this problem, the corresponding solution of the 
modulation system must depend on the self-similar variable r = £/£ alone, which reduces 
the modulation system to the ordinary differential equations 

dv ' 

K-r)-^ = 0, 2 = 1,2,3. (79) 
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Figure 5: Left: Riemann invariants behaviour in the similarity modulation solution for the 
flat-bottom zero-friction case ; Right: corresponding undular bore profile ?/(£)• 



The boundary conditions for ( J791 follow from the matching conditions (J73l) . ( 1741) using the 
initial condition ( J7H1) : 



r 
r 



r 



?"2 
?'2 



r 3 



r 3 



-A 



0. 



30) 



where r 1 * 1 are self-similar coordinates (speeds) of the leading and trailing edges, £ = r 1 ^. 
Taking into account the inequality r\ < r 2 < r3 one obtains the well-known modulation 
solution of Gurevich and Pitaevskii (1974) (see also Fornberg and Whitham 1978) in the 
form 

r\ = —A , r 3 = , r 2 = -mA , (81) 

£ /a . /~,\ f\Ar/-i x 2m(l — m)K(m) 
j = u 2 (-A, -mA, 0) = 2A[(1 + 



(82) 



E(m) — (1 — m)K(m) 

This modulation solution (I8~T|) . (152]) (see Fig. 5a) represents the replacement, due to averag- 
ing over the oscillations, of the unphysical formal three-valued solution of the dispersionless 
KdV equation (i.e. of the Hopf equation) which would have taken place in the absence of 
the dispersive regularisation by the undular bore. We see that fl82l describes an expansion 
fan in the characteristic (£, C)-plane and thus is a global solution. Substituting flHTj) . fl82|) 
into the travelling wave solution ( |70l) one obtains the asymptotic wave form of the undular 
bore (see Fig. 5b), which then can be readily represented in terms of the original physical 
variables using the relationships ( 1681) . 

The equations of the trailing and leading edges of the undular bore are determined from 
(|82|) by putting m = and m = 1 respectively 



Y = r ~ = ^(-A, 0, 0) = -6A , ^ = r+ = v 2 (-A, -A, 0) = 4A . 



(83) 



The leading solitary wave amplitude is r/ = 2(r 3 — r x ) = 2A, which is exactly twice the height 
of the initial jump. This corresponds to the amplitude of the surface elevation a = 3h A (see 
(|68|) ). Note that, to get the leading solitary wave of the same initial amplitude ao as for the 
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separate solitary wave considered in Section 4, one should use the jump value A = a /3h , 
which of course is just 2 A, where A = 3/ioA/2 is the initial discontinuity in the surface 
elevation. 



6.3 Structure of the undular bore front 

We are especially interested in the behaviour of the modulation solution (IHTj) . (1821) in the 
vicinity of the leading edge £ = £ + (C)- This behaviour is essentially determined by the 
manner in which the pair of characteristics corresponding to the velocities vi and vi merge 
into a multiple eigenvalue v + of the modulation system at £ = £ + (C)- 

First, one can readily infer from the modulation solution (1HTT) . ( 1821) that the phase velocity 
c = — 2(ri + r 2 + r 3 ) = 2A(1 + m) > t> 2 (— A, — m, 0) for m < 1 and c = v 2 for m = 1. Thus, 
any individual wave crest generated at the trailing edge of the undular bore moves towards 
the leading edge, i.e. for any crest m — > 1 as ( — > oo. Thus, for any particular wave crest, 
except for the very first one, the solitary wave 'status' is achieved only asymptotically as 
C -> oo. 

Without loss of generality we assume in this section that A = 1 in ( IHTi) . ( 1821) . First, as 
we have already mentioned, the characteristic family r 2 : d£/d( = t> 2 is an expansion fan in 
the £, ( - plane, 

T 2 : £ = C 2 C, (84) 

parameterised by a constant C 2 , —6 < C 2 < 4 . Next, in fl82|) we make an asymptotic 
expansion of f 2 (— 1, — m, 0) for small (1 — m) <C 1, to get 

2(1 -m)ln(16/(l - m)) ~ r+ (85) 

or, with logarithmic accuracy, 

(r + -«/ C )«l: l-^ 2 , n[1/(T+ ^ /c)r (86) 

Next, expanding t>i(— 1, — m, 0) for (1 — m) <^ 1 and using fl86l) we get the asymptotic 
equation for the characteristics family Ti, 

^ = Vl = r + + ( r + _ £/ C ) + 0(1 _ m ) , (87) 
which is readily integrated to leading order to give 

Ti: ^r+C-y, (88) 

where C% > is an arbitrary constant 'labeling' the characteristics; C\ = corresponds to 
the leading edge of the undular bore. This asymptotic formula (IHHj) is valid as long as ( 3> 1. 
The behaviour of the characteristics belonging to the families Ti and T 2 near the leading 
edge is shown in Fig. 6a. 

Next, expanding the equation for the third characteristic family, T 3 : d£/d( = v 3 (— 1, — m, 0) 
for (1 - m) < 1, we get on using flHBI) 

|- M1 A T "-;/C)) +0 < T+ -^- (89 > 
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Figure 6: Characteristics behaviour for the similarity modulation solution near the leading 
edge £ + (C) : ( a ) families IV d£/d( = v\ and r 2 : £ = C2C, (b) family T 3 : d£/d( = t> 3 . 

Integrating fl89|) we obtain to first order 

T 3 : f ~ C 3 - 9(C) , (90) 

where A0 = f S (C 3 /r+) = 0, (91) 

C*3 being an arbitrary constant. The asymptotic formula f)90p is valid as long as g{C)/C^ <C 
1. Since the characteristics 1^ intersect the leading edge £ = r + £ we must indicate their 
behaviour outside the undular bore. It follows from the matching condition (|74|) and the 
limiting structure fl34l) of the characteristic velocities of the Whitham system, that the 
characteristics from the family 1^ match with the Hopf equation characteristics d£/d( = 6r 
carrying the value of the Riemann invariant r = corresponding to still water upstream the 
undular bore. Therefore, the sought external characteristics are simply vertical lines £ = C3. 
The qualitative behaviour of the characteristics from the family 1^3 is shown in Fig. 6b. 

It is clear from the asymptotic behaviour of the characteristics that the edge characteristic 
£ = t + ( corresponding to the motion of the leading solitary wave intersects only with 
characteristics of the family T 3 carrying the Riemann invariant value r 3 = into the undular 
bore domain. Since, according to the matching condition (1801) . r 3 = everywhere along the 
edge characteristic one can infer that the leading solitary wave motion is completely specified 
by its amplitude at ( = 0. Indeed, in this case, the leading edge represents a genuine multiple 
characteristic of the modulation system, along which the Riemann invariant r + = r 2 = r\ is 
a constant. Given the constant value of r\ = —1 for the solution (|82p . one infers that the 
amplitude of the lead soliton of the self-similar undular bore, r] Q = 2(r 3 — r + ) = 2 is also a 
constant value. Thus, in the undular bore evolving from an initial jump, the leading solitary 
wave represents an independent soliton of the KdV equation. Of course, this fact follows 
directly from the modulation solution ([8*21 but now we have established its meaning in the 
context of the characteristics, which will play an important role below. 

Next we discuss the structure of the undular bore front in the case when the initial 
profile //(£, 0) is not a simple jump discontinuity, and instead has the form of a monotonically 
decreasing function, for instance, (— O 1 ^ 2 when £ < and r)(£, 0) = for £ > 0. In that case, 
the modulation solution for the undular bore no longer possesses x/t-similarity as in the 



23 




Figure 7: a) Leading edge £ + (C) of non-self-similar undular bore as an envelope of pairwise 
merging characteristics from the families d£/d( = v% and d^/dC, = v<i] b) behaviour of the 
Riemann invariants in non-self-similar modulation solution with r 3 = 0. 

jump resolution case and, as a result, the speed (and therefore, the amplitude) of the lead 
solitary wave is not constant. For instance, for the afore-mentioned square-root initial profile 
the amplitude of the lead solitary wave grows as ( 2 (see Gurevich, Krylov and Mazur 1989, 
or Kamchatnov 2000). Clearly, such an amplitude variation is impossible if the leading edge 
£ + (C) was a regular characteristic carrying a constant value of the Riemann invariant r + . As 
discussed above, however, the GP matching conditions f!73|) -f J77|) admit another possibility; 
the leading edge curve is the envelope of the characteristic families r^: d£/d( = V\ and T 2 : 
d^/dQ = t> 2 merging when m = 1. This configuration is shown in Fig. 7a. In this case, the 
behaviour of the modulus m in the vicinity of the leading edge is given by the asymptotic 
formula found in Gurevich & Pitaevskii (1974): 

where the function r + (Q ^ constant is assumed to be known. Another specific feature of 
this (general) configuration is that dri^/d^ — > ±oo as £ —>■ £ + (see Fig. 7b - also found in 
Gurevich & Pitaevskii 1974, see also Kamchatnov 2000), which is in drastic contrast with 
similarity solution (see Fig. 6a). This particular difference was discussed in relation with 
undular bores in the KdV-Burgers equation in Gurevich and Pitaevskii (1987). 

In summary, we see from (1921) that the structure of the modulation solution in the vicin- 
ity of the leading edge of an undular bore defined as a characteristic envelope is qualitatively 
different compared to that for the similarity case (see (J85l) ). The more general (but qual- 
itatively similar to (192jl ) asymptotic formula which takes into account small perturbations 
due to a variable topography and bottom friction will be derived later. At the moment, 
it is important for us that in this configuration, when the leading edge is a characteristic 
envelope rather than just a characteristic, the value r + , and thus, the leading solitary wave 
amplitude are allowed to vary. 

The analysis of the corresponding modulation solution in Gurevich, Krylov and Mazur 
(1989) showed that, while in the case of an initial jump the wave crests generated at the 
trailing edge reach the leading edge (and therefore, transform into solitary waves) only 
asymptotically as t — > oo, for the more general case of decreasing initial data each wave 
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crest generated at the trailing edge reaches the leading edge in finite time and replaces 
(overtakes) the existing leading solitary wave. This process is manifested as a continuous 
amplitude growth of the (apparent) leading solitary wave. As in classical soliton theory, 
an alternative explanation of the leading solitary wave amplitude growth can be made in 
terms of the momentum exchange between the "instantaneous" leading solitary wave and 
solitary waves of greater amplitude coming from the left. Indeed, as the rigorous analysis of 
Lax, Levermore and Venakides showed (see Lax, Levermore and Venakides (1994) and the 
references therein), the whole modulated structure of the undular bore can be asymptotically 
described in terms of the interactions of a large number of KdV solitons initially 'packed' 
into a non-oscillating large-scale initial profile. 

This latter interpretation is especially instructive for our purposes. Our point is that 
the specific cause of the enhanced soliton interactions resulting in amplitude growth at the 
leading edge is not essential; it can be large-scale spatial variations of the initial profile as 
just described, but it could also equally well be an effect of small perturbations in the KdV 
equation itself. Indeed, in the weakly perturbed KdV equation, the local wave structure 
of the undular bore must be described to leading order by the periodic solution (ITU]) of 
the unperturbed KdV equation, so if one assumes the GP boundary conditions analogous to 
(I75|) - ([77]) for the perturbed modulation system ([27]) . one invariably will have to deal with 
one of the two possible types of the characteristics behaviour (shown in Figs. 7a and 8a) 
in the vicinity of the leading edge of the undular bore, because this qualitative behaviour 
is determined only by the structure of the GP boundary conditions and by the associated 
asymptotic structure of the characteristic velocities of the Whitham system for (1 — m) <C 1, 
which are the same for both unperturbed and perturbed modulation systems. Next, we will 
show that, by using the knowledge of this qualitative behaviour of the characteristics, one 
is able to construct the asymptotic modulation solution for the undular bore front in the 
presence of variable topography and bottom friction even if the full solution of the perturbed 
modulation system is not available. 

6.4 Gurevich-Pitaevskii problem for perturbed modulation sys- 



We investigate now how the GP matching problem applies to the perturbed modulation 
system (]27j) . As in the original GP problem, we postulate the natural physical requirement 
that the mean value (U) is continuous across the undular bore edges, which represent free 
boundaries and are defined by the conditions m = (trailing edge X = X~(T)) and m — 1 
(leading edge X = X + (T)). Also, we consider propagation of the undular bore into still 
water, hence (U)\ x =x+(t) — 0. Now, using the explicit expression (fTSl) for (U) in terms of 
complete elliptic integrals and calculating its limits asm->0 and m — > 1 one has 



tern 



X = X~(T) : 
X = X + (T) : 



A2 — A; 

A 2 = A 



1 • 



(U) 
(U) 



-Ai = u, 
-A 3 = 0, 



(93) 



where u(X,T) is solution of the dispersionless perturbed KdV equation ([7]), i.e. 



u T + 6uu x = F(T)u - G(T)u 2 , 



(94) 
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with the boundary conditions 



U \ 6# / hdT ^) = 2T~ A ° if T < T °' M ( T ' 6g / = ° if T > T ° ' ^ 

where r = —x /^gh . The boundary conditions fl95|) correspond to a discontinuous initial 
surface elevation A(x,t) at x = — x , obtained by using transformations ([3]) and ([6]) where 
one sets t = 0. As earlier, A = ao/(3h ) is the value of the discontinuity in A, chosen in 
such a way that the amplitude of the lead solitary wave in the undular bore was exactly ao 
in the flat-bottom zero-friction region (see Section 6.2). 

This free-boundary matching problem is then complemented by the kinematic conditions 
explicitly defining the boundaries X = X ± (T). These are formulated using the multiple 
characteristic directions of the perturbed modulation system (|27|) in the limits as m — > 
and m -> 1 (cf. fl7S} - (177jl ). 

<1X ~ V~(X-,T), d -^r = V + {X\T), (96) 



where = ^(w, A - , A - ) = vs(u, A~, A - ), (97) 

V + = v 2 (\ + , A+, 0) = ^(A+, A+, 0) , (98) 

and \~ = \ 2 (X-,T) =\ 3 (X~,T), A + = X 2 (X + ,T) = X 1 (X + ,T). (99) 

Thus, for the perturbed KdV equation the leading and trailing edges of the undular bore are 
defined mathematically in the same way as for the unperturbed one, albeit for a different 
set of variables. 



6.5 Deformation of the undular bore front due to variable topog- 
raphy and bottom friction 

Finally we study the effects of gradual slope and bottom friction on the leading front of the 
self-similar expanding undular bore described in Sections 6.2, 6.3. The result will essentially 
depend on the relative values of the small parameters appearing in the problem. We note 
that in general there are three distinct relevant small parameters, 

e = — < 1 , 5 = max(/i,) < 1, C D < 1 (100) 

The first small parameter is determined by the ratio of the equilibrium depth in the flat 
bottom region, to the distance from the beginning of the slope region to the location of the 
initial jump discontinuity in the surface displacement. This measures the typical relative 
spatial variations of the modulation parameters in the undular bore when it reaches the 
beginning of the slope. The second and third parameters are contained in the KdV equation 
(OQ) itself and measure the values of the slope and bottom friction respectively. In terms of 
the transformed variables appearing in ([7]), |-F(T)| ~ 5, \G(T)\ ~ Co (see OH])). Generally 
we assume 5 ~ Co (the possible orderings 5 Co or Co "C 5 can be then considered as 
particular cases). 
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To obtain a quantitative description of the vicinity of the leading edge of the undular 
bore we perform an expansion of the Whitham modulation system (|2"T|) for (1 — m) <C 1. 
We first introduce the substitutions 



X i {X,T)=\+{T) + l i {X,T), Vt = V + + v[, p l =p + + p' l , i 



1,2. 



where X = X~ 



X 



-4A + , p* 



A ^9 

-F(T)A+ + -G(T)(A+) 2 . 



(101) 
(102) 



Since A2 > Ai, v 2 > v\ one always has l 2 > h, v 2 > v[. Assuming X / X + Cl<^>l - m < 1 
and using that A3 = to leading order in the vicinity of the leading edge (see the matching 
condition (|93|) ). we have from asymptotic expansions of (1281) - (I3T]) as (1 — m) <C 1 



v[ = M l {l 2 -l 1 ) = -2 
v' 2 = M 2 (l 2 - h) = -2 



ln(16/(l -m)) 



1 - 



l + |(l-m) ln(16/(l -m))_ 

ln(16/(l-m)) 
1 -i(l-m) ln(16/(l-m)) 



(ia-ii), 
(/2-/1), 



(103) 



1 + ln 



4 
15 



p' 2 = N 2 {l 2 -h) 



2A + ln 



5 + ln 



-16A+ 



-16A+ 

h-h 
-16A+ 



3A 



G (I 2 -Ii 



(104) 



4 

+ 15 



Gki 2 -ii). 



Naturally, v[ and p£ vanish when l 2 Now, substituting ( 11011) . ( 11021) into the modulation 
system (1271) we obtain 



9/- 

cff 1 dX dT ( V+ + ^ = P + +P- 



1,2. 



On using the kinematic condition fl96l) at the leading edge, this reduces to 

= 1,2. 



d\ + ,dk + , / 



(105) 



(106) 



There are two qualitatively different cases to consider: 

(i) lim^o \dk/dX\ < 00, i = 1, 2 (Fig. 8a) 

(ii) lim^o Itfli/dX"! = 00, i = 1, 2 (Fig. 8b) 

The case (i) implies that to leading order (I106j) reduces to 

d\+ 



dT 



(107) 



which, together with the kinematic condition dX + /dT = — 4A + , defines the leading edge 
curve X + (T). One can observe that this system coincides with (T43|) . (142]) defining the 
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x + x x- x 

a) b) 



Figure 8: Riemann variables behaviour in the vicinity of the leading edge of the undular 
bore propagating over gradual slope with bottom friction (a) Adiabatic variations of the 
similarity GP regime, 5 <C e, Cd <C e; (b) General case, 5 ~ Cd ~ 6. 



motion of a separate solitary wave over a gradual slope with bottom friction. Its integral 
expressed in terms of original physical x, t- variables is given by fj49|) . Therefore, in the case 
(i) the lead solitary wave in the undular bore to leading order is not restrained by interactions 
with the remaining part of the bore and behaves as a separate solitary wave. Physically this 
case corresponds to adiabatic deformation of the similarity modulation solution f[8"Tj) . f[52"j) 
and implies the following small parameter ordering : 5 <C e, Cd *C e. 

Next, we study the structure of this weakly perturbed similarity modulation solution in 
the vicinity of the leading edge. The next leading order of the system (11061) yields 



dX 



1,2, 



(108) 



that is 



dh _Ni dl 2 _N 2 
~dX~~^~i dX~~M~ 2 
Subtraction of one equation (11091) from another with account of the relationship l 2 — l\ 
— A + (l — m) leads consistently to leading order to the differential equation for 1 — m 



(109) 



8(1 — m) 
dX 



F(T) 16G(T) 



15 



In 



16 



m 



-i 



This equation should be solved with the initial condition 

l_ m = at X = 0. 
Elementary integration gives with the accuracy 0(1 — m) (cf. (155]) ) 



[1 — m) In ■ 



16 



1 — m 



X + -X 
-A+ ' 



(110) 



[111) 



(112) 



This formula determines the dependence of the modulus monT and X (as long as 1 — m 
I)- 
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Now, we make use of the solution A + of equation (11071) given by (1471) with C = 
4/(3ga h Q ) (see (]4"5]) ). Under supposition that the integral f x h~ 3 dx diverges as h — > 0, 
so that the turbulent bottom friction plays an essential role in the undular bore front be- 
haviour (see Section 4 for a similar approximation for an isolated solitary wave), we obtain 
for h <C ho 

At last, if the bottom topography is approximated by the dependence (152"]) . we get with the 
same accuracy 



16 64 
[l-m)\n—^ = \-C D 
1 — m 15 



2 + 



(3a — 1)5 \h 



(X+-X), (114) 



where a > 1/3. The second term in square brackets tends to zero as h — > 0. However, the 
region where it can be neglected may be very narrow because of smallness of the parameter 
5. We recall that in this formula X + is given by ( 149]) and X is defined by ([3]) in terms of 
the original physical independent variables x and t. 

Summarising, if the conditions 8, Cp <C e are satisfied, the lead solitary wave of the 
undular bore behaves as an individual (noninteracting) solitary wave adiabatically varying 
under small perturbation due to variable topography and bottom friction. The modulation 
solution in the vicinity of the leading edge also varies adiabatically, however, its qualitative 
structure considered in Section 6.4 (see Figs 5,6) remains unchanged. 

In a sharp contrast with the described case of adiabatic deformation of an undular bore 
front is case (ii) when the second term in the left-hand side of fll06j) contributes to the 
leading order, i.e. to the motion of the leading edge itself. Namely, we have 

d\~^~ i , dlrj 

—. = p+ + v '—L i = l,2. 115 
dT H l dX 

Now dX + 1 dT 7^ p + which means that the amplitude of the lead solitary wave a = — 2A + varies 
essentially differently compared to the case of an isolated solitary wave. Indeed, the term 
p + in the right-hand side of ( 11151) is responsible for local adiabatic variations of the solitary 
wave while the term v[dli/dX describes nonlocal parts of the variations associated with 
the wave interactions within the undular bore. Using asymptotic formulae (11031) implying 
v' 2 > 0, v[ < 0, and the condition \im^^ \dli^/dX\ = oo along with li > l±, it is not 
difficult to show that this nonlocal term is always nonnegative , i.e. the lead solitary wave 
in the undular bore propagating over a gradual slope with bottom friction always moves 
faster (and, therefore, has greater amplitude) than an isolated solitary wave of the same 
initial amplitude in the beginning of the slope. Indeed, as we have shown in Section 5, the 
presence of a slope and bottom friction always result in "squeezing" the cnoidal wave, hence 
increasing momentum exchange between solitary waves in the vicinity of the leading edge 
of the undular bore and acceleration of the lead solitary wave itself. The situation here is 
qualitatively analogous to that described in Section 6.4 where the general global modulation 
solution for the unperturbed KdV equation was discussed. Similarly to that case, the leading 
edge now represents a characteristic envelope - a caustic (otherwise we are back in the case 
(i) implying d\ + /dT = p + ) (see Fig. 6a). 
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Unlike the case of adiabatic variations of the leading edge, determination of the function 
A + (T) requires now knowledge of the full solution of the perturbed modulation system (127)) 
with the matching conditions (I93p . While the analytic methods to construct such a solution 
for inhomogeneous quasilinear systems are not available presently, it is instructive to assume 
that d\ + /dT — p + is a known function of T and to study the structure of the solution in 
close vicinity of the leading edge. With an account of the explicit form (11031) of the velocity 
corrections, equations ( 11151) assume the form 



dX 

dk 
dX ' 



d\+/dT - p+ \ 1 1 .„ 

2{l 2 -h) [ln[16/(l - m)] ^ 4 l } 



dX + /dT-p^ 



ln[16/(l -m)] 4 



+ -(1 - m) 



2(/ 2 - h) 

Taking the difference of H 1 1 6 1) and fll 17[) we transform it to the form 
9(1 -m) d\+/dT-p+ 1 



dX (A+) 2 (1 -m)ln[16/(l - m)\ 

This equation can be readily integrated with the initial condition (11 1 in to give 



(116) 
(117) 

(118) 



:i-m) 



In 



16 1 
1 - m + 2 



2(dX + /dT - p H 
(A + ) 2 



-(X+-X). 



;ii9) 



This solution coincides with the asymptotic formula ()92l) for the behaviour of the modulus 
in the vicinity of the leading edge of the undular bore in general unperturbed GP problem 
[T6] but instead of the derivative d\ + /dT in ( 192)) we have the difference d\ + / dT — p + (which 
is always positive as we have established). 



7 Conclusions 

We have studied the effects of a gradual slope and turbulent (Chezy) bottom friction on the 
propagation of solitary waves, nonlinear periodic waves and undular bores in shallow-water 
flows in the framework of the variable-coefficient perturbed KdV equation. The analysis has 
been performed in the most general setting provided by the associated Whitham equations 
describing slow modulations of a periodic travelling wave due to the slope, bottom friction 
and spatial nonuniformity of initial data. This modulation theory, developed in general form 
for perturbed integrable equations in Kamchatnov (2004) was applied here to the perturbed 
KdV equation and allowed us to take into account slow variations of all three parameters 
in the cnoidal wave solution. The particular time-independent solutions of the perturbed 
modulation equations were shown to be consistent with the adiabatically varying solutions 
for a single solitary wave and for a periodic wave propagating over a slope without bottom 
friction obtained in Ostrovsky & Pelinovsky (1970, 1975) and Miles (1979, 1983a). It was 
shown, however, that the assumption of zero mean elevation used in these papers for the 
description of slow variations of a cnoidal wave, ceases to be valid in the case when the 
turbulent bottom friction is present. In this case, a more general solution was obtained 
numerically improving the results of Miles (1983b). 
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Further, the derived full time-dependent modulation system was used for the descrip- 
tion of the effects of variable topography and bottom friction on the propagation of undular 
bores, in particular on the variations of the undular bore front representing a system of 
weakly interacting solitary waves. By the analysis of the characteristics of the Whitham 
system in the vicinity of the leading edge of the undular bore, two possible configurations 
have been identified depending on whether the leading edge of the undular bore represents a 
regular characteristic of the modulation system or its singular characteristic, i.e. a caustic. 
The first case was shown to correspond to adiabatically slow deformations of the classi- 
cal Gurevich-Pitaevskii modulation solution and is realised when the perturbations due to 
variable topography and bottom friction are small compared with the existing spatial non- 
uniformity of modulations in the undular bore (which is supposed to be formed outside the 
region of variable topography/bottom friction). In the case when modulations due to the 
external perturbations are comparable in magnitude with the existing modulations in the 
undular bore, the leading edge becomes a caustic, and this situation was shown to corre- 
spond to enhanced solitary wave interactions within the undular bore front. These enhanced 
interactions have been shown to lead to a "nonlocal" leading solitary wave amplitude growth, 
which cannot be predicted in the frame of the traditional local adiabatic approach to prop- 
agation of an isolated solitary wave in a variable environment. As we mentioned in the 
Introduction, one of our original motivations for this study was the possibility to model a 
shoreward propagating tsunami as an undular bore. In this context, we would suggest that 
the second scenario described above is the more relevant, which has the implication that 
the growth, and eventual breaking of the leading waves in a tsunami wavetrain, cannot be 
modeled as a local effect for that particular wave, but is determined instead by the whole 
structure of the wavetrain. 
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Appendix A: Derivation of the perturbed modulation system 

We express the integrand function in the right-hand side of (j21J) in terms of the ^-variable 
(HI: 

(2A, - Sl - U)R = 8G/i 3 - [8G\i + 4(F + 2 Sl G)]/i 2 

+ [4(F + 2s 1 G)X i + 2s 1 (s 1 G + F)}fi-2 Sl ( Sl G + F)\i. ^ ' 

Then we obtain with the use of f[T3"j) . (|14|) . and ffTol) the following expressions: 
(fi) = y (f fid9 — \ (f n^-dn = \ ^ 1 ' 



LJ LJ d\i L J 2y/-P(fj.) Lds 2 
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where J is a known integral 



A 3 



a/ (A3 - - A 2 )(/x - Ai) d\i 



(122) 



15 



(A 3 - A!) 5/2 [(l -m + m 2 )E(m) - (1 - m)(l - m/2)K(m)}, 



K[m) and E(m) being the complete elliptic integrals of the first and second kind, respec- 
tively. The derivatives of / with respect to A« are also known table integrals (Gradshtein & 
Ryzhik 1980): 



dl 
3A7 



dl 

dX, 



dl 

ox: 



A3 (A3 _^ )( ^_ A2 ) 



A2 



jJL - Aj 



^V / A 3 -Ai[(A 2 + A 3 - 2A X )E - 2(A 2 - Xi)K], 



a 3 

A 2 



(A 3 - //)(// - Ai) 



/i - A 2 



-i\/A3-Ai[(A 3 - Xi)K + (Ai + A 3 - 2A 2 )£], 



(123) 



A3 



A 2 



(/iA 2 )(/z - A 



A 3 - n 



—djj, 



i v / A 3 -A 1 [(2A 3 - Ai - A 2 )£ - (A 2 - X^K). 



We can easily express the Sj-derivatives in terms of Aj derivatives by differentiation of the 
formulae (see (fIBl) ) 



s i — Ai + A 2 + A3, s 2 — AiA 2 + A1A3 + A 2 A3, S3 — AiA 2 A3 
and solving the linear system for differentials. Simple calculation gives 

dXt f-l) 3 - fc 



(124) 



(125) 



ds k n^^Ai-Aj)' 

Then, combining f!123j) and (I125p . we obtain the derivatives dl/dsi and hence the expressions 
2 



~L ~ 15 
1 dl 1 



(A3 - Xx) 



(si - 3s 2 )§ - h\ 2 - Ax)(A 2 + A 3 - 2Ax) 

A Z 



E 



L dsi 6 
1 dl 1 
Lds^ " ~2 



2 Sl — + s 1 X 1 + X(- A 2 A 3 



(126) 



(A 3 -AO— + Ai 
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To complete the calculation of the right-hand side of (|24|) . we need also expressions 

1 2(A 2 -A 1 )f, 



dL/d\ 1 y z E 

L 2(A 3 - X 2 )(l-m)K 

dL/d\ 2 ~ E-(l-m)K ' 

L _ 2(A 3 - X 2 )K 
dL/dX 3 ~ E-K ' 



(127) 



Collecting all contributions into perturbations terms, we obtain the Whitham equations in 
the form 

f)\ f)\ 

df + v ^ = P* = C>[F(T)A t - G{T)Bl (128) 
where Cj, A,, Bj and v j, j = 1, 2, 3 are specified by formulae (1281) - (130]) . 
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